Odd triplet superconductivity in superconductor/ferromagnet multilayered structures. 
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We demonstrate that in multilayered superconductor- 
ferromagnet structures a non-coUinear alignment of the mag- 
netizations of different ferromagnetic layers generates a triplet 
superconducting condensate which is odd in frequency. This 
triplet condensate coexists in the superconductors with the 
conventional singlet one but decays very slowly in the ferro- 
magnet, which should lead to a large Josephson effect between 
the superconductors separated by the ferromagnet. Depend- 
ing on the mutual direction of the ferromagnetic moments the 
Josephson coupling can be both of and n type. 

It is well known that ferromagnetism and supercon- 
ductivity are antagonistic phenomena; ferromagnetism 
strongly suppresses superconductivity. This suppression 
in superconductor/ferromagnet (S/F) layered structures 
is caused mainly by the exchange interaction. This means 
in particular that the singlet Cooper pairs cannot pene- 
trate into the F layers over a noticeable length, since the 
exchange energy is by orders of magnitude larger than 
the Cooper binding energy. Thus, the singlet pairs are 
destroyed by the exchange field because the spins of the 
electrons cannot be antiparallel anymore. This suppres- 
sion of superconductivity can be reduced in S/F struc- 
tures if the ferromagnetic layers are ultrathin, but, gen- 
erally, the destruction of the singlet superconductivity by 
a homogeneous exchange field can hardly be avoided. 

The situation may be different if the spins of the super- 
conducting pairs are parallel to each other. It is clear that 
such a triplet superconductivity is not sensitive to the ex- 
change field and the coexistence of the superconductiv- 
ity and ferromagnetism becomes possible. Unfortunately, 
the triplet pairing is a rather exotic phenomenon and has 
been observed until now only in superfluid He^ and in a 
superconducting material Sr2Ru04 [1]. It is expected to 
be very sensitive to disorder [1,2], which makes its obser- 
vation even more difficult. In order to satisfy fermionic 
commutation relations the condensate function, even in 
frequency, should be odd in the momentum of the pair 
and this is the reason why it is so sensitive to impuri- 
ties. Another possibility was suggested by Berezinskii 
long ago [3]. He conjectured that the triplet supercon- 
ductivity might be possible if the condensate function 
were even in the momentum but odd in the frequency. 
Attempts to find conditions for the existence of such an 
odd superconductivity were done in several works much 
later [4] but the results were not encouraging (in Ref. [4] 
a singlet pairing odd in frequency and in the momentum 
was considered). 



In this Letter, we demonstrate that the odd triplet 
superconductivity is not exotic at all and is inavoidable 
in multilayered S/F structures with conventional super- 
conductors if the directions of the magnetization of the 
different F-layers are neither parallel nor antiparallel to 
each other. In this case the triplet condensate (TC) can 
easily penetrate into the ferromagnetic layers over long 
distances and result in supercurrents through the ferro- 
magnets. The conditions for the realization of the odd 
triplet superconductivity do not seem to be problematic 
from the experimental point of view and we hope that 
proper measurements will be done in the nearest future. 

A generation of a triplet condensate by a non- 
homogeneous magnetization (domain walls) and its pen- 
etration into a ferromagnet has been discussed recently 
[5,6]. In the structures considered in these works only 
changes in the conductivity of the system were analyzed. 
However, the detection of the triplet component in such 
structures is a quite difficult task. In contrast, we pre- 
dict in this Letter a new type of superconductivity (odd 
triplet) in S/F structures and discuss how to identify it 
experimentally. To be specific, we consider a S/F mul- 
tilayered structure (Fig. 1) in which the new supercon- 
ducting state might be observed. In this state, the su- 
perconductivity in the S regions is caused by the singlet 
component (SC) and takes place in the plane of the lay- 
ers, while the transverse superconductivity through the 
F layers is mainly due to the TC. Moreover, the relation 
between the condensate current Is and a phase difference 
if depends in a crucial way on chirality of the magnetic 
moment M in space. 

The multilayered S/F structure we consider is repre- 
sented in Fig.l. The (in-plane) magnetizations M of the 
neighboring F layers are not parallel to each other and 
the angle between them is 2a. In order to achieve such 
a non-coUinear alignment one can employ an exchange- 
biased spin- valve or ferromagnets with strong anisotropy 
and different easy-axis of magnetization. It will be shown 
that in such a structure the TC may arise if the thick- 
ness of the superconducting layers 2ds is less or compa- 
rable with the coherence length — \/Ds/2ttT. The 
TC penetrates into the F layer over a long distance 
= \JF)f/2i^T and ensures a Josephson coupling be- 
tween the nearest S layers. In the case under consid- 
eration the relation between the condensate current Is 
and the phase difference Lp has the conventional form 
Is — Ic sin (/9. However, the sign of the critical current 
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Ic depends on the chirality, namely, it is positive if the 
rotation angle 2aj of the magnetization at the Si layer 
has the same sign for neighboring S layers (S, and St+i) 
and it is negative if the signs of the rotation are opposite 
for neighboring S layers (ajaj+i < 0). 
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FIG. 1. Schematic diagram of a S/F multilayered struc- 
ture. On top: the magnetization M in different F layers. The 

M vector drawn by solid lines correspond to chiralitics of the 
same sign. Opposite chirality is obtained if, for example, the 
magnetization of the right F layer is rotated in the opposite 
direction (dashed line). 

Now, let us explain how these results were obtained. 
We consider the simplest case of a dirty system for which 
the Usadel equations can be applied to calculate the con- 
densate function /, i.e., we assume that the condition 
A < J < is fulfilled, where A is the order parameter 
in S (in the ferromagnet A = 0), J is the exchange energy 
and T is the momentum relaxation time. However, even 
in this case the problem remains too complicated and we 
need to make further assumptions. First, we assume that 
the thickness of the F layers 2dF exceeds the length ^t- 
In this case the Josephson coupling is determined by an 
overlap of exponentially decaying condensate functions 
induced by the nearest S layers. We also assiimc that 
the temperature is close to the critical temperature T* 
of the structure (|r - r*|/T* < 1). In this case the 
Usadel equation can be linearized and what remains to 
do is to find the condensate function / {x) by solving 
the linearized Usadel equation in the S and F regions. 
The Usadel equation for the matrix condensate function 
/ takes the form [7,8] 

DsdlJ-2\uj\f = 2iA, Slayer (1) 



Dpdl^f - 2|w|/ + i Jsgnw { [cts, /]+ cos a+ 
+T3[o'2,/]-sina} = 0, Flayer (2) 

Here Ds^f are the diffusion coefficients in the S and F 
layers respectively, A = zf2 (EiiraA, uj = ■7rT(2n-|- 1) is the 
Matsubara frequency, J is the exchange energy (strictly 
speaking, J is a vector directed along the magnetization 
vector M), fj and o-j are the Pauli matrices in the Nambu 



and spin space, respectively. The brackets [...]± denote 
an anti-commutator and a commutator. We assume that 
the M vector lies in the {y, z) plane and a is the an- 
gle between the z— axis and M. Eqs.(l-2) have to be 
supplemented by boundary conditions. In the case of a 
perfect S/F interface (the reflection coefficient is very 
small) they have the form [9] 



f{ds + 0) ~ f{ds - 0) = 0, 
ldJ{ds + 0)-dJ{ds-Q) = Q, 



(3) 



where 7 = ap/as, (yp.s are the conductivities in the F 
and S regions (for simplicity we do not take into account 
a dependence of the conductivity ap on the spin direc- 
tions). If the SF interface resistance Rb is finite the r.h.s 
of the first equation equals Jb^jdxfids + 0), where 7(, is 
proportional to Rb [10]. For simplicity we set 7^ = 0. 
A generalization of the results for the case of a finite jb 
does not lead to qualitative changes. 

The condensate matrix function / has following form 



/(a;)=?(/i(a;)fi + /2(x)f2) 



(4) 



where /i^2 are matrices in the spin space. In the F layers 
they can be written as 



/l = Bi (x) CTi + B2 (x) (72, 

h = Bo (x) + B3 {x) as . 



(5) 



In the S layers the functions /i^2 have the same form, but 
the coefficients B should be replaced by A. The functions 
Bi (x) have the form: Bi{x) ^ bi cxp(=pK(a; =p ds)). Sub- 
stituting these expressions for Bi{x) into Eqs. (1-2) one 
obtains a set of algebraic equations that determine the 
eigenvalues k. In the S region all condensate components 
(SC and TC) arc decoupled and there is only one solu- 
tion for k: — 2\uj\/Ds- In the F region the SC and 
TC arc coupled by the exchange interaction (if a is not 
zero) and there are three eigenvalues for k. Two of them 
are equal to k± = kj{1 ± i), and determine a fast decay 
of the condensate in the F layers ( dashed line in Fig. 
2), here kj = \J J j Dp . This result demonstrates the 
well known short-range penetration of the superconduct- 
ing condensate into the ferromagnet (see for example the 
review article [11]). However, the third solution for k 
is completely different and is given by = y^lwI/Dp. 
Thus, we see that the huge exchange energy J entering 
the solution for k± is replaced by a small energy to of the 
order of T. At large distances it is the solution with 
that determines the penetration of the superconducting 
condensate into the F regions (see Fig. 2). It corresponds 
to the triplet component and, as it has been discussed in 
the introduction, there is no wonder that the exchange 
field does not influence it. 

Having determined the solutions for k, we can write the 
solutions for Bi (x). For example, the coefficient Bi{x) 
may be written at x > in the form 
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-^i(^) = ^^P [-i^ii^ - ds)] , (6) 

i 

where bn = bi^^ , bi± and Kj = Ki^^ , Ki± . The coefRcients 
bi± are related to the coefficients 63^ that determine the 
magnitude of the SC by bi± = ^sinasgno; • 63^. 

What remains to be done is to determine the coefR- 
cients Oj, bi using the boundary conditions (3) for an ar- 
bitrary angle a. In a general case the expressions for the 
coefficients are cumbersome, and for clarity we present 
here only results for small a. In zero order in a only co- 
efficients in the function /2 differ from zero. For example, 
for 63^ we obtain 



2i\Lu\ 1 + ^k±/{ks tanh^^s) 

where ks = ■\/2\lo\/Ds, Os = Ksds- This solution corre- 
sponds to those obtained earlier (see [11] and references 
therein) where the case of a parallel or antiparallel mag- 
netization alignment was analyzed. It is valid in the lim- 
iting cases ^5 << 1 and >> 1 ( if T is close to the 
critical temperature T*). 

If the magnetization vectors in the neighboring F lay- 
ers are inclined with respect to each other by an angle 
2a, the situation changes qualitatively. First, the func- 
tion /i that describes the TC is no longer zero because 
the terms proportional to bi± in the expression for Bi{x) 
(Eq. (6)) are related to &3±, and hence they are finite in 
the main approximation. Secondly, in order to fulfill the 
matching conditions (3), one has to take into account 
the first long-range term in Eq.(6). Using the boundary 
conditions, we obtain for the coefficient bi,^ the following 
expression 

A Kj tanh 9 s 

— sm a 2 ; o 

w cosh 6s\'yK+/Ks + tanh^sl (kw tanh^g + Ks/7) 

(8) 

Contrary to the SC determined by the solution -63(0;) 
and ^3(2;), the TC is an odd function of x. This is seen 
from Eqs. (6), (8) and the fact the a has different sign 
at a; > ^5 and x < —ds- According to the boundary 
conditions (3) (continuity of the condensate fmiction and 
current conservation), the TC induced in the F layer pen- 
etrates into the superconductor. The corresponding so- 
lution in the S layer has the form: Ai{x) = ai sinh(K5a;). 
In the limit T « J, the TC penetrates into the F layer 
over the length of the order 1/kt = \/Dp/27tT which is 
much greater than the SC penetration length 1/k±. One 
can see from Eq. (8) that the amplitude bia, of the long- 
range TC is an odd function of the Matsubara frequency 
oj and is symmetric in the momentum space as in the 
case of the TC which arises in a ferromagnet with a non- 
homogeneous M near a S/F interface [5]. This new type 
of the condensate, odd in w and even in the momenta 



p, has been proposed by Bcrczinksii [3] in order to ex- 
plain the pairing mechanism in He^ (it was proven later 
that the condensate in He^ is in fact even in w and odd 
in p). The solution wc present here corresponds to this 
hypothetical pairing, which means that we have found a 
concrete realization of the idea. It follows from the ge- 
ometry of the structure we consider that the odd triplet 
superconductivity we found exists in the transverse di- 
rection. 

It is seen from Eqs. (7) and (8) that the amplitudes 
of the SC and TC at the S/F interface are comparable 
if << and a and 6s are of the order of 1. If 
the latter condition is not satisfied (that is, the thickness 
of the S layer 2ds is large in comparison with ^s), the 
TC decays exponentially in S and its amplitude is small. 
Therefore we calculate the Josephson current between the 
neighboring S {Si and 5*2) layers in the case 0s <C 1 and 
assuming that the condition « < 2dF is fulfilled. 
Then, the Josephson current is due to the overlap of the 
TC induced near each S/F interface. In this case the TC 
in F is described by the expression 

ftrip{x) = in (g) (Ti6ia, exp(-Kt^(x - ds))+ 
-I- 5 • if 1 (g) (Ti • S+biuj exp(K„(a; - ds - 2^^)) (9) 

Here, the matrix of a gauge transformation S = 
cos{ip/2) + if3sin((p/2) allows us to take into account 
the phase difference tp between the neighboring S lay- 
ers (we assume that the phase of the 5*1 layer is zero). 
The coefficients have different signs (61^^ = — ^i^^) if the 
magnetization M at both S layers rotates in the same 
direction, and bi^i = bi^i in the opposite case (different 
signs of chiralities). The condensate current Is between 
and S2 is given by the formula 

Is = {LyL,)cjF{-KiT/Ae)Tr{h ® ^o) fdj (10) 

From Eqs. (9-10), we get Is = IcSintp, where 
Ic = -{2'kT / eRF)2dF ^ Kw^Iw^Iw exp(-2rfF«a>) (H) 

OJ 

and Rf = 2dF/{{I'yI'z)<^F) is the resistance of the 
F layer in the normal state. Substituting Eq.(8) into 
Eq.(ll), we find for the critical current in the limit -C 1 

_ AT /A\' [K.,ds)-(dt KT)v--'''-^'^' smm ■ smn2 
^ TT \^^Sods + li^j{l + i)lnsonnTds + lhY 

(12) 

where k,so = \/2ttT/Ds- 

If the magnetization vector at each S layer rotates in 
the same direction (ai = a^), then the critical Josephson 
current Ic is positive (0-contact). If M rotates at Si and 
S2 in opposite directions (ai = —a2), then the critical 
current Ic is negative (tt— contact). In the latter case a 
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phase difference tt is established between neighboring S 
layers in a multilayered S/F structure. We would like 
to note that the mechanism of the tt— contact considered 
here is completely different from that suggested in Ref. 
[12] and observed in Ref. [13]. In our case, the negative 
critical current is caused by the TC, but not by the SC 
as in Ref. [12], and, in addition, is realized only if the 
chiralities at the Si and 5*2 are different. The possibility 
of switching between the and 7r-states by changing the 
angle a may find applications in superconducting devices. 

The effect described above exists for all temperatures 
below the critical temperature T* which is determined 
from the self-consistency equation. In the main approxi- 
mation in the small angle a it agrees with that obtained 
in Refs. [10], [11] for the case of the parallel orientation 
of the magnetization. A correction to T* due to a small 
misalignment of the magnetizations in neighboring F lay- 
ers, i.e. due to the TC, is proportional to sin^ a, but we 
are not interested in this correction here. We note that 
the critical temperature T* for the case of arbitrary a 
was analyzed in Ref. [14]. However, the form of the con- 
densate function presented in that work is not correct 
because the authors started from an equation different 
from Eq. (2 ) (instead of the commutator, they wrote 
the product Ja ■ f ). As a result, the long range triplet 
component was completely lost. 
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FIG. 2. The normalized absolute value of the SC (dashed 
line) in the S and the F layers. The solid line repre- 
sents the absolute value of the TC in the F layer. Note 
that the SC (TC) is an even (odd) function of x. Here 
ds/^T ~ d.s/£,s ~ 0.3, ds/Cj ~ 5, "f — 0.1 and sina — 0.5. 
We have assumed a finite SF interface resistance and set 
76 = 0.5. 



superconductors separated by the ferromagnet is possible 
and the critical current can be measured. The Josephson 
contact can be both of and 7r-type depending on the 
arrangement of the magnetic moments. We hope that 
the effects considered in this paper will be observed ex- 
perimentally in the nearest future. 
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In conclusion, we have predicted a new type of 
superconductivity. It was demonstrated that in 
superconductor-ferromagnet structures a non-coUinear 
alignment of the exchange fields in the ferromagnetic lay- 
ers generates a triplet component of the superconducting 
condensate and this component is odd in the frequency. 
The odd triplet condensate penetrates into the ferromag- 
net over long distances and is not sensitive to impurities. 
In the structure we suggest, a Josephson effect between 
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